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Modified  embedded  atom  method  study  of  the  mechanical  properties  of  carbon  nanotube 

reinforced  nickel  composites 

Jamal  Uddin,1,2’*  M.  I.  Baskes,3  S.  G.  Srinivasan,2,4  Thomas  R.  Cundari,1,2  and  Angela  K.  Wilson1,2 
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2 Center  for  Advanced  Scientific  Computing  &  Modeling,  University  of  North  Texas,  Denton,  Texas  76203,  USA 
3 Materials  Science  and  Technology  Division,  Los  Alamos  National  Laboratory,  MS  G755,  Los  Alamos,  New  Mexico  87545,  USA 
4 Department  of  Materials  Science  &  Engineering,  University  of  North  Texas,  Denton,  Texas  76203,  USA 
(Received  12  October  2009;  published  11  March  2010) 

We  report  an  atomistic  simulation  study  of  the  behavior  of  nanocomposite  materials  that  are  formed  by 
incorporating  single-walled  carbon  nanotubes  (SWCNTs),  with  three  different  diameters,  and  a  multiwalled 
carbon  nanotube  (MWCNT)  into  a  single-crystal  nickel  matrix.  The  interactions  between  carbon  and  nickel 
atoms  are  described  by  a  modified  embedded  atom  method  potential.  Mechanical  properties  of  these  nanocom¬ 
posite  materials  are  predicted  by  atomistic  calculations  and  compared  with  that  of  fee  nickel  and  pristine  CNTs. 

Our  simulations  predict  that  all  Ni/CNT  composites  studied  in  this  work  are  mechanically  stable.  Their  elastic 
properties  depend  on  the  volume  fraction  and  diameter  of  embedded  CNTs.  The  single-crystal  Young’s  modu¬ 
lus  (Elf)  of  Ni/SWCNT  composites  exhibit  a  large  increase  in  the  direction  of  CNTs  alignment  compared  to 
that  of  a  single-crystal  nickel.  However,  a  moderate  but  gradual  decrease  is  seen  for  E22  and  £33  in  the 
transverse  directions  with  increase  in  CNT  diameters.  As  a  consequence,  Ni/SWCNTs  show  a  gradual  decrease 
for  the  poly  crystalline  Young’s,  bulk  and  shear  moduli  with  the  increasing  CNT  diameters  and  volume  frac¬ 
tions.  These  reductions,  although  moderate,  suggest  that  enhancement  of  mechanical  properties  for  polycrys¬ 
talline  Ni/SWCNT  nanocomposites  are  not  achievable  at  any  CNT  volume  fraction.  The  Ni/MWCNT  com¬ 
posite  with  high  CNT  volume  fraction  shows  the  highest  increase  in  En.  Unlike  the  E22  and  E33  for  Ni/ 

SWCNTs,  there  is  a  significant  increase  in  the  E22  and  the  E33  for  Ni/MWCNT.  As  a  result,  poly  crystalline 
Ni/MWCNT  composites  show  slight  increase  in  the  elastic  properties.  This  suggests  that  nickel  nanocompos¬ 
ites  with  enhanced  mechanical  properties  can  be  fabricated  using  large  volume  fractions  of  larger  diameter 
MWCNTs.  Depending  on  type,  alignment  and  volume  fraction,  Ni/CNT  composites  show  varying  degrees  of 
elastic  anisotropy  and  Poisson’s  ratio  compared  to  pure  Ni.  Simulation  predicts  strong  adhesion  at  the  Ni/CNT 
interface  and  a  significant  interfacial  stress  transfer  between  CNT  and  Ni  matrix. 

DOI:  10.11 03/Phy sRe vB . 8 1 . 1 04 1 03  PACS  number(s):  62.23.Pq,  81.40.Jj,  68.35.Gy,  62.20.de 


I.  INTRODUCTION 

Carbon  nanotubes  (CNTs)  are  among  the  stiffest  and 
strongest  materials.1  Their  extremely  small  size,  high  aspect 
ratios,  and  extraordinary  structural  and  chemical  stability 
make  CNTs  a  leading  candidate  for  reinforcement  fillers  in 
nanocomposite  materials.2-6  As  a  result,  there  has  been  an 
increasing  scientific  and  technological  interest  in  developing 
nanocomposites  by  dispersing  CNTs  into  a  metal  matrix  or 
by  laser  depositing  metal  matrix  around  the  CNTs.  Such 
composites  combine  low  density  with  high  strength  and  stiff¬ 
ness,  and  enhanced  electrical  and  thermal  properties.  These 
are  particularly  attractive  for  aerospace  and  other  transporta¬ 
tion  applications2-4  where  there  is  a  renewed  emphasis  on 
increasing  fuel  efficiency  by  “light  weighting”  the  transport¬ 
ers. 

CNT-metallic  composites  offer  distinct  advantages  over 
polymeric  composites  for  high-strength  applications.5,6  Ku- 
zumaki  et  al.  were  the  first  to  study  a  aluminum/CNT 
composite.7  Preparation  and  application  of  magnesium/CNT 
composites  was  also  reported.8,9  The  mechanical  strength  of 
nickel/CNTs  nanocomposites  was  studied  recently  by  Sun  et 
al. 10  Some  progress  has  also  been  made  in  the  application  of 
atomistic  simulations  to  predict  the  elastic  and  the  mechani¬ 
cal  properties  of  such  composite  materials.  For  example, 
molecular-dynamics  (MD)  simulations  have  been  success¬ 


fully  applied  to  predict  the  elastic  properties  of  polyethylene/ 
CNT  composites11  and  other  polymeric  CNT  composites.12 
However,  despite  their  potential  superiority  in  mechanical 
properties,  the  development  of  CNT-metal  matrix  composites 
remains  in  its  infancy.  Development  of  metal  matrix  compos¬ 
ites  will  be  promoted  from  atomistic  insights  on  the  response 
of  their  mechanical  properties  to  chemical  and  material 
modifications. 

Nickel  is  the  primary  ingredient  of  numerous  aerospace 
superalloys,13,14  and  therefore,  is  an  attractive  target  for  fab¬ 
rication  of  ultrastrength  materials  with  CNTs.  This  research 
investigates  the  mechanical  properties  of  pristine  CNTs,  a  fee 
nickel  matrix,  and  Ni/CNT  composites  at  different  CNT  vol¬ 
ume  concentrations  using  quasistatic  atomistic  simulations. 
Elastic  properties  of  solids  are  important  because  they  relate 
to  various  fundamental  solid-state  properties.  They  deter¬ 
mine  the  response  of  the  crystal  to  external  forces,  as  char¬ 
acterized  by  elastic  constants,  bulk  modulus,  shear  modulus, 
Young’s  modulus,  and  Poisson’s  ratio.  Therefore,  elastic 
properties  play  an  important  part  in  determining  the  strength 
of  the  materials.  We  have  calculated  these  properties  and  the 
elastic  anisotropy  of  Ni/CNT  nanocomposite  materials.  Inter¬ 
atomic  interactions  are  described  using  a  modified  embedded 
atom  method  (MEAM)  initially  developed  by  Daw  and 
Baskes15,16  and  later  modified  by  Baskes.17 
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TABLE  I.  Parameters  for  the  MEAM  potentials  of  Ni,  C,  and  Ni-C.  The  parameters  are  the  cohesive 
energy  Ec  (eV),  the  equilibrium  nearest-neighbor  distance  re  (A),  the  exponential  decay  factor  for  the  uni¬ 
versal  energy  function  a ,  the  scaling  factor  for  the  embedding  energy  A,  the  four  exponential  decay  factors 
for  the  atomic  densities  f£l\  the  four  weighting  factors  for  the  atomic  densities  t^\  and  the  density  scaling 
factor  p°.  Angular  screening  parameters  for  the  MEAM  potentials  for  binary  systems. 


MEAM  potentials 

Angular  screening  parameters 

Ni 

C 

Ni-C 

c 

^max 

c 

^min 

Ec 

4.45 

7.37 

4.82 

Ni-Ni-Ni 

2.8 

0.8 

re 

2.49 

1.54 

2.01 

Ni-C-Ni 

2.8 

2.0 

a 

4.99 

4.38 

4.82 

C-Ni-C 

2.8 

2.0 

A 

1.10 

1.49 

Ni-Ni-C 

2.8 

2.0 

2.45 

4.26 

Ni-C-C 

2.8 

2.0 

/3d) 

1.50 

5.00 

C-C-C 

2.8 

2.0 

£(2) 

6.00 

3.20 

/30) 

1.50 

3.98 

t(0) 

1.00 

1.00 

5.79 

7.50 

1.60 

1.04 

3.70 

-1.01 

P° 

1.00 

1.00 

II.  COMPUTATIONAL  DETAILS 


A.  Theoretical  formalism 


In  MEAM,  the  total  energy  of  a  system  is  given  in  the 
following  form: 


F ’i(Pi)  +  ^  E  SijVijiRij) 
K*i) 


where  Fi  is  the  embedding  function  for  an  atom  i  embedded 
in  a  background  electron  density  pp  SLj  is  the  screening  func¬ 
tion,  and  <Pij(Rij)  is  the  pair  interaction  between  atoms  i  and 
j  separated  by  a  distance  Rir  For  the  calculations  of  energy, 
the  functional  forms  of  Ft  and  cp^  should  be  given.  The  back¬ 
ground  electron  density  at  each  atomic  site  is  computed  by 
combining  several  partial  electron-density  terms  for  various 
angular  contributions  with  weight  factors  t^h\h  =  0-3)  (direc¬ 
tionality  of  bonding).  Each  partial  electron  density  is  a  func¬ 
tion  of  atomic  configuration  and  atomic  electron  density.  The 
atomic  electron  densities  pa^h\h  =  0-3)  are  given  as 

pa(h>(R')  =  p()  exp[-  0h)(R/re  -  1 )], 

where  p0,  the  atomic  electron-density  scaling  factor,  and  0h\ 
the  decay  lengths,  are  adjustable  parameters,  and  re  is  the 
nearest-neighbor  (NN)  distance  in  the  equilibrium  reference 
structure.  A  specific  form  is  given  to  the  embedding  function 
F i  but  not  to  the  pair  interaction  cp^.  Instead,  a  reference 
structure  where  individual  atoms  occupy  perfect-crystal  atom 
sites  is  defined  and  the  potential  energy  per  atom  of  the  ref¬ 
erence  structure  is  estimated  from  the  zero-temperature  uni¬ 
versal  equation  of  state  by  Rose  et  al.ls  The  value  of  the  pair 
interaction  is  evaluated  from  the  known  values  of  the  poten¬ 
tial  energy  per  atom  and  the  embedding  energy,  which  is  a 
function  of  the  nearest-neighbor  distance. 


In  the  original  MEAM,17  only  the  first  nearest-neighbor 
interactions  were  considered.  The  second  (2NN)  and  more 
distant  neighbor  interactions  were  neglected  by  using  strong 
many-body  screening  function.  The  second-nearest-neighbor 
interactions  have  been  included  in  the  modified  formalism  by 
adjusting  screening  parameters  Cmin  to  decrease  the  many- 
body  screening  effect.  In  addition,  a  radial  cutoff  function  is 
also  applied  to  reduce  calculation  time.  The  MEAM  for  an 
alloy  system  is  based  on  the  MEAM  potentials  for  the  com¬ 
ponent  elements.  The  2NN  MEAM  formalism19  gives  14  in¬ 
dependent  model  parameters  for  pure  elements:  four 
(Ec,re,B,d)  for  the  universal  equation  of  state,  seven 

for  the  electron  density,  one 
(A)  for  the  embedding  function,  and  two  (Cmin, Cmax)  for  the 
many-body  screening.  The  details  of  the  MEAM  formalism 
have  been  published  in  the  literature. 17,20,21 

B.  Simulation  methods 

The  MEAM  potential  parameters  that  describe  the  inter¬ 
atomic  interactions  for  Ni,  C,  and  Ni-C  composites  are  listed 
in  the  Table  I.  It  treats  nickel,  carbon,  and  the  Ni-C  interac¬ 
tions,  and  allows  the  formation  of  chemical  bonds  with  ap¬ 
propriate  atomic  hybridization.  Similar  MEAM  potential  pa¬ 
rameters  have  been  shown  to  model  the  catalytic  growth  of 
CNTs  on  nickel  nanocatalysts  well.22  We  performed  molecu¬ 
lar  static  calculations  at  T=0  K.  Periodic  boundary  condi¬ 
tions  were  applied  along  all  three  directions  of  the  supercells. 
The  Ni  matrices  and  Ni/CNT  composites  were  relaxed  to 
minimize  their  energies. 

The  Ni/CNT  composites  were  simulated  using  three 
single- walled  nanotubes  (SWCNTs)  with  varying  diameters 
and  a  multiwalled  (MWCNT)  nanotube.  We  chose  CNTs 
with  zigzag  structures  for  this  study.  SWCNTs  with  (5,0), 
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®  O 

m  m 

SWCNT(5,0)  SWCNT(10,0) 

(a)  (b) 

FIG.  1.  (Color  online)  Space  filling  unit  cells  of  (n,0)  CNTs. 
SWCNT  and  MWCNT  refer  to  single- walled  and  multi  walled  CNT. 
(a)  SWCNT(5,0),  (b)  SWCNT(10,0),  and  (c)  SWCNT(15,0). 
MWCNT  in  (d)  is  built  by  assembling  SWCNT(5,0), 
SWCNT(10,0),  and  SWCNT(15,0)  nanotubes  coaxially.  The  ap¬ 
proximate  tubular  height  for  the  unit  cell  is  4.26  A  for  all 
nanotubes. 

(10,0),  and  (15,0)  chiral  indices  and  tubular  diameters  of 
3.97  A,  7.86  A,  and  11.77  A,  respectively,  were  studied. 
These  CNTs  were  embedded  in  a  supercell  containing  37  044 
nickel  atoms.  The  atomic  structure  and  the  crystallographic 
orientation  of  the  SWCNTs  and  supercells  are  shown  in  Figs. 
1  and  2.  This  supercell  has  21X21X21  unit  cells  and  is 
approximately  7.4  nmX7.4  nmX7.4  nm  in  three  perpen¬ 
dicular  directions. 

As  the  CNT  volume  fraction  is  an  important  variable  in 
determining  the  composite  mechanical  properties,  we  used 
two  different  concentrations  of  each  CNTs  for  Ni/CNT  com¬ 
posites.  The  three  supercells  had  a  CNT  volume  fraction  of 
0.83%,  1.85%,  and  3.37%  when  single  CNTs  were  used. 
When  four  CNTs  are  embedded  in  each  supercell,  the  CNT 


id) 


FIG.  2.  (Color  online)  Atomistic  structures  and  the  crystallo¬ 
graphic  orientations  of  the  supercells,  (a)  Fully  minimized  pure 
nickel  (face-centered-cubic  lattice),  (b)  nickel  matrix  with  an  em¬ 
bedded  MWCNT,  (c)  nickel  matrix  with  four  embedded  MWCNT, 
(d)  full  extent  of  a  MWCNT  embedded  in  the  nickel  supercell. 
Periodic  boundary  conditions  are  applied  for  pure  nickel  and  Ni / 
CNT  matrix  along  the  three  dimensions. 


volume  fractions  were  increased  to  3.33%,  7.40%,  and 
13.46%  for  (5,0),  (10,0),  and  (15,0)  SWCNTs,  respectively. 
Since  the  Ni  matrix  does  not  penetrate  the  hollow  inside  of 
CNT,  we  treated  the  CNTs  as  solid  beams  for  volume  frac¬ 
tions  calculations. 

The  MWCNTs  were  built  by  inserting  SWCNTs  into  each 
other  coaxially  while  maintaining  their  individual  periodic¬ 
ity.  For  composites  made  with  MWCNTs,  the  CNT  volume 
fractions  were  similar  to  the  volume  fractions  of 
SWCNT(15,0).  For  example,  Fig.  2(d)  shows  MWCNT  that 
fits  the  periodicity  of  fee  Ni  matrix.  The  wall  thickness  of 
SWCNTs  were  assumed  to  be  0.32  nm  based  on  electron- 
density  distribution  calculations23  and  similar  methods.24,25 

To  simulate  the  Ni/CNTs  nanocomposites  with  low  CNT 
volume  fractions,  four  different  single-crystal  Ni  supercells 
containing  a  hollow  cylindrical  void  at  the  center  were  built. 
For  studying  nanocomposites  with  high  CNT  volume  frac¬ 
tions,  four  different  single-crystal  Ni  supercells  were  built 
with  four  cylindrical  voids  approximately  equidistant  to  each 
other.  The  cylindrical  void  runs  from  one  face  to  the  opposite 
face  of  the  supercell  with  the  cylinder  axis  aligned  parallel  to 
the  supercell  [001]  direction  ( a  axis).  The  radii  of  these  voids 
were  chosen  to  fit  SWCNT(5,0),  SWCNT(10,0), 

SWCNT(15,0),  and  the  MWCNT,  leaving  an  appropriate 
amount  of  space  for  Ni-C  interatomic  bonding.  The  nano¬ 
composites  were  formed  by  embedding  the  CNTs  into  the 
voided  Ni  supercells.  The  resulting  composites  contain  a  ho¬ 
mogeneous  distribution  of  aligned  CNTs.  Two  simulation 
cells  with  embedded  CNTs  are  shown  in  Figs.  2(b)  and  2(c). 

Mechanical  properties  and  stress-strain  relationships  were 
determined  from  the  relaxed  configurations.  Elastic  constants 
for  Ni  and  Ni/CNTs  were  calculated  from  strain-energy- 
applied  strain  relationships  by  fitting  the  strain  energy  per 
volume  to  a  parabola.  Strains  of  0.1%  deformation  around 
the  fully  minimized  structures  were  chosen  to  maintain  de¬ 
formations  within  the  linear  elastic  regime.  For  each  defor¬ 
mation  step  (0.0002  step  size),  we  minimized  system  energy 
by  relaxing  positions  of  all  atoms.  Nine-independent  elastic 
constants  (c11?  C2 2,  C33?  Ci2>  £*23’  c3b  c44>  c55?  and  £*66)  and 
three  independent  constants  (cn,  c12,  and  c44)  of  Ni/CNTs 
and  pure  fee  Ni  were  determined  in  this  study.  The  direction- 
dependent  Young’s  modulus  ( E ),26  which  measures  the  stiff¬ 
ness  of  a  material,  was  calculated  from  the  slope  of  the 
strain-stress  curves  for  pristine  CNTs.  Young’s  modulus  for 
CNTs  with  only  longitudinal  response  were  calculated  with 
0.0002  increments  in  the  applied  strain.  For  Ni  matrix  and 
Ni/CNT  systems,  E  was  calculated  from  the  appropriate 
single-crystal  elastic  constants. 

III.  RESULTS  AND  DISCUSSION 
A.  Pristine  CNTs 

Table  II  shows  that  the  Young’s  moduli  of  individual 
CNTs  obtained  in  our  work  agree  well  with  numerous  pub¬ 
lished  experimental28,30,31,33  and  theoretical  values.22,24,27,43 
Most  of  the  experimental  measurements  of  E  were  carried 
out  for  MWCNTs.  Also,  Young’s  moduli  of  SWCNTs  show  a 
large  variation  from  0.78  to  3.6  TPa  and  that  only  a  few 
experimental  results  are  available.33,35,37,44,45  Our  calculated 
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TABLE  II.  Young’s  modulus  (£n)  (in  GPa)  for  single-walled  (SWCNTs)  and  a  multiwalled  (MWCNT)  carbon  nanotubes  obtained  from 
our  ME  AM  Ni-C  potential  are  tabulated  and  compared  with  the  literature  values. 


This  work 

Previous  calculations 

Experiments 

820  (Xiao  et  al,  MD)a 

900  (Agrawal  et  al. ,  MD)b 

270-900  (Yu  et  al.f 

5500  (Yakobson,  MD)d 

790  (Li  et  al.f 

552  [SWCNT(5,0)] 

974  (Lu,  MD)f 

900  (Demcyzek  et  al.f 

832  [SWCNT(10,0)] 

1325  (Cheng  and  Gao,  MD)h 

816  ±410  (Salvetat  et  al.)' 

892  [SWCNT(15,0)] 

1238  (Jin  and  Yuan,  Forcc-Ficld)J 

1200  (Tombler  et  al.f 

740  (MWCNT) 

1220  (Hernandez  et  al.,  TB-MD)1 

4700  (Tu,  LDA)n 

764  (Zhou  et  al.,  LCAO-MO  cluster)0 

965,  972  (Ogata  and  Shibutani,  TB,  DFT)P 

1470  (Pullen,  GGA^ 

1060  (Li  and  Chou,  MM)r 

899(5,0),  916(10,0),  923(15,0)  (Meo  and 

Rossi,  MM-FEM)8 

1250  (Krishnan  et  al.)m 

Reference  22. 
bReference  27. 
cReferenee  28. 
dReference  29. 
Reference  30. 


Reference  24. 
gReference  31. 
hReference  32. 
Reference  33. 
^Reference  34. 


kReference  35. 
Reference  36. 
“Reference  37. 
nReference  38. 
°Reference  39. 


pReference  40. 
qReference  41. 
rReference  32. 
Reference  42. 


E  of  832  GPa  for  SWCNT(10,0)  agrees  with  the  value  cal¬ 
culated  by  Xiao  et  al  (820  GPa),  Agrawal  et  al.11  (900  GPa), 
and  Lu  (975  GPa).24  Our  results  from  the  same  CNT  is  also 
in  reasonable  agreement  with  the  values  of  965  and  972  GPa 
calculated  using  tight-binding  and  density-functional  theory, 
respectively,  by  Ogata  and  Shibutani.40  On  the  contrary, 
Hernandez  and  coworkers36  reported  a  much  greater  value  of 
1220  GPa  for  this  nanotube  using  tight-binding  MD.  These 
large  differences  may  be  attributed  to  differences  in  the  CNT 
volume  estimation  due  to  variations  in  effective  CNT  wall 
thickness.  In  spite  of  extensive  studies,  the  variations  in  wall 
thickness  of  CNTs  remain  an  unresolved  question  46 

B.  Ni/CNT  composites: 

1.  Single-crystal  elastic  constants 

A  necessary  condition  for  mechanical  stability  of  a  struc¬ 
ture  is  that  its  strain  energy  must  be  positive  against  arbitrary 
but  small  homogeneous  elastic  deformations.  For  an  ortho¬ 
rhombic  crystal,  this  imposes  the  following  constraints:47 
c11  +  c22>2c12,  C22  +  C33>2c23,  C1 1  +  C33  >  2c31,  Cu>  0  (/ 

—  1 — 6),  and  Cjj+c 22~t~c33~t~2c^2~i~2c23~i~2c3i^:>0.  These  crite¬ 
ria  simplify  to  cn>0,  c44>0  c11  +  2c12>0,  and  cn-cu 
>0  for  cubic  crystals  such  as  fee  Ni.  Hence,  we  examined 
the  mechanical  stability  of  the  Ni/CNT  composites  using  ap¬ 
propriate  single-crystal  elastic  constants  for  composites  with 
both  low  and  high  CNT  volume  fractions.  The  calculated 
values  for  elastic  constants  c ^  are  summarized  in  Table  III. 

For  single-crystal  Ni,  the  MEAM  calculated  values  for 
cn,  c12,  and  c44  agree  well  with  experiments  48,49  The  shear 
elastic  constants  cn-c12  and  cn  +  2c12  determine  stability. 


Physically,  c44  plays  a  subtle  role  and  represents  the  resis¬ 
tance  of  the  system  to  a  shear  deformation  such  as  in  the 
(100)  plane  in  the  [010]  direction.  Similarly,  (cn-c12)/2 
represents  the  resistance  to  shear  deformation  across,  for  ex¬ 
ample,  a  (110)  plane  in  the  [110]  direction.  Similar  observa¬ 
tions  hold  for  an  orthorhombic  crystal.  From  the  Table  III 
data,  it  is  observed  that  for  Ni/CNT  composites,  the  cu  ( i 
=  1-6),  cn+c22-2c12,  c22+c33-2c23,  cu+c33-2c13,  and 
cn+c22+c33  +  2c12+2c23  +  2c13  are  all  positive,  and  thus  it 
can  be  reasonably  concluded  that  Ni/CNTs  composites  are 
mechanically  stable.  This  is  an  encouraging  result  for  the 
fabrication  of  nanocomposite  materials  by  embedding  the 
CNTs  in  metal  matrices. 

For  Ni/CNTs,  cn  show  moderate  increase  in  composites 
with  low  CNT  volume  fractions  and  significantly  large  in¬ 
crease  in  composites  with  high  CNT  volume  fractions  com¬ 
pared  to  that  of  pure  Ni.  These  increases  depend  on  the  type 
and  diameter  size  of  the  embedded  nanotubes.  For  example, 
composites  Ni/SWCNT(5,0),  Ni/SWCNT(10,0),  and  Ni/ 
SWCNT(15,0)  show  an  increase  of  5,  10,  and  12  GPa,  re¬ 
spectively,  for  low  CNT  volume  fraction  and  19,  39,  and  42 
GPa,  respectively,  for  high  CNT  volume  fraction.  Ni/ 
MWCNT  shows  the  largest  increase  in  cn  than  for  any  of  the 
Ni/SWCNTs.  In  this  case,  cn  shows  an  increase  of  30  GPa 
and  109  GPa,  respectively,  for  low  and  high  CNT  volume 
fractions  compared  to  that  of  Ni.  In  the  transverse  directions, 
c22  or  c33  show  gradual  decrease  for  Ni/SWCNTs  compared 
to  pure  Ni.  However,  in  these  cases,  composites  with  higher 
CNT  volume  fractions  show  a  larger  decrease  compared  to 
composites  with  low  CNT  volume  fractions.  Interestingly, 
the  Ni/MWCNT  nanocomposite  shows  very  small  or  no 
changes  for  c21  or  c33. 


104103-4 

4 


MODIFIED  EMBEDDED  ATOM  METHOD  STUDY  OF  THE... 


PHYSICAL  REVIEW  B  81,  104103  (2010) 


TABLE  III.  Single-crystal  elastic  constants  (in  GPa),  and  derived  quantities  cn  +  c22-2c12,  c22+c33-2c23,  ^11  +  c33-2c13,  and  cn 
+  C22+c33  +  2c12+2c23  +  2c13  for  pure  fee  nickel  and  the  Ni/SWCNTs  and  Ni/MWCNT  composites  with  low  and  high  CNT  volume  fractions. 


Ni(fcc) 

Ni/SWCNT(5,0) 

Ni/SWCNT(10,0) 

Ni/SWCNT(15,0) 

Ni/MWCNT 

Low  High 

Low  High 

Low  High 

Low  High 

cn 

c22 

c33 

C12 

c23 

C13 

c44 

c55 

c66 

c11  +  c22~2c12 
c22  +  c33~^c23 
cll  +  c33~^cl3 

Cll  +  C22  +  C33  +  ^C12  +  ^C23  +  ^C13 


247.1 

251.9 

266.1 

247.1 

246.0 

243.6 

247.1 

246.0 

243.9 

147.5 

146.6 

144.0 

147.5 

148.3 

150.5 

147.5 

148.0 

144.9 

125.4 

123.7 

118.6 

125.4 

121.4 

112.3 

125.4 

121.4 

112.3 

199.3 

204.7 

221.7 

199.3 

201.3 

220.1 

199.3 

195.5 

186.5 

1626.2 

1630.3 

1632.5 

256.7 

284.0 

258.8 

245.4 

241.5 

243.2 

245.2 

240.6 

243.4 

145.5 

139.7 

142.7 

147.9 

148.8 

147.7 

145.5 

140.1 

142.7 

118.5 

99.6 

114.5 

121.9 

113.1 

120.7 

121.9 

113.1 

120.7 

211.1 

246.1 

216.6 

210.9 

244.4 

216.8 

194.8 

184.5 

191.2 

1625.1 

1623.3 

1611.6 

289.4 

276.9 

356.1 

233.4 

245.7 

243.9 

233.6 

246.3 

247.5 

129.7 

142.9 

132.8 

147.5 

149.1 

149.6 

129.5 

143.6 

133.8 

83.8 

118.2 

109.5 

109.1 

118.4 

104.7 

109.1 

118.4 

104.8 

263.8 

236.8 

334.3 

263.6 

236.1 

294.9 

172.0 

193.8 

184.5 

1569.8 

1640.0 

1728.7 

In  the  composites  with  SWCNTs,  c12  and  c13  decrease  for 
both  low  and  high  CNT  concentrations,  which  depend  on  the 
CNT  diameters  and  the  CNT  volume  fractions.  A  larger  de¬ 
crease  is  seen  for  composites  with  larger  CNT  diameters.  For 
example,  in  the  case  of  Ni/SWCNT(15,0)  with  high  CNT 
volume  fraction,  the  value  for  c12  and  c13  is  found  to  be  18 
GPa  less  than  that  of  Ni.  The  c23  on  the  other  hand,  remain 
similar  to  that  of  Ni,  with  a  few  exceptions,  where  a  minute 
increase  is  seen.  For  the  composites  with  the  MWCNT,  c12 
and  c13  decrease  for  both  low  and  high  CNT  volume  frac¬ 
tions.  For  c23,  however,  a  small  increase  is  seen  for  both 
volume  fractions.  The  values  of  c44,  c55,  and  c66  decrease  in 
the  composites  compared  to  pure  Ni  (125  GPa)  and  depend 
on  the  CNT  concentration  in  the  composites.  A  similar  trend 
is  seen  for  c55  and  c66,  although  these  show  less  variation 
than  c44.  The  Ni/MWCNT  composites  also  show  similar 
CNT  concentration  dependence  for  these  elastic  parameters, 
however,  the  decrease  in  c44  is  less  in  a  relative  sense.  Note 
that  c44,  c55 ,  and  c66  are  equivalent  to  the  shear  modulus  (G). 

To  study  the  angular  anisotropy  of  the  elastic  constants  in 
Ni/CNT,  we  performed  calculations  in  which  the  CNT  is 
rotated  about  its  tubular  axis.  Elastic  constants  from  these 
calculations  are  similar  to  the  original  unrotated  systems. 
While  lack  of  anisotropy  merits  further  investigation,  the  oc¬ 
currence  of  a  regular  interfacial  structural  pattern  seems  un¬ 
likely  because  this  interface  forms  between  the  large  periph¬ 
eral  surface  of  the  CNT  that  runs  through  the  entire  7.4  nm 
length  of  the  unit  cell  and  the  irregular  curvature  of  the  Ni 
matrix  in  the  perpendicular  directions. 

The  Young’s  modulus  (E)  along  the  normal  direction  of 
the  lattice  plane  is  defined  as  the  ratio  of  the  stress  to  the 
strain  along  that  direction.  En  (E001),  E22  (Emo),  and  E33 
(E0iq)  were  calculated  for  cubic  and  orthorhombic  lattices 
along  the  principal  ( a ,  b ,  and  c)  axes.  Poisson’s  ratio  v(n,m) 
of  an  elastic  solid  for  any  two  specified  general  orthogonal 
unit  vectors  n  and  m  is  defined  as  -s(m)/ s(n),  where  s(m)  is 
the  transverse  strain  in  the  direction  m  and  s(n )  is  the  axial 


strain  in  the  direction  n  due  to  a  uniaxial  tension  of  the 
material  along  the  direction  n.  With  the  exception  of  isotro¬ 
pic  elastic  materials,  Poisson’s  ratio  for  an  anisotropic  elastic 
material  depends  on  the  choice  of  n  and  m  and  was  calcu¬ 
lated  from  the  elastic  compliance  values. 

Figure  3  presents  the  calculated  values  of  ^n*  e22 ,  and 
E33  for  Ni  and  Ni/CNT  composites  for  two  different  CNT 
volume  fractions.  For  pure  Ni,  we  calculated  E  of  136.9  GPa 
along  the  major  directions.  This  agrees  reasonably  well  with 
an  experimental  value  of  120  GPa.49  For  composites,  in  the 
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FIG.  3.  (Color  online)  Calculated  Young’s  modulus  (GPa)  for 
fee  Ni  and  Ni/CNT  composites  with  two  different  CNT  fractions. 
The  F001  (En)  is  the  Young’s  modulus  for  uniaxial  loading  parallel 
to  CNT  axis,  the  F010  (E22)  or  F100  (F33)  are  the  Young’s  moduli  in 
the  directions  transverse  to  the  CNT  axis  in  single-crystal  systems, 
and  Fpolycrystalline  is  the  Young’s  modulus  for  poly  crystalline  aggre¬ 
gates.  Composites  with  low  CNT  fractions  are  designated  as  CL  and 
composites  with  high  CNT  fractions  are  designated  as  CL. 
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FIG.  4.  (Color  online)  Calculated  Poisson’s  ratio  (v)  for  fee  Ni 
and  Ni/CNT  composites  with  low  and  high  CNT  volume  fractions. 

longitudinal  direction  (the  direction  of  CNT  alignments),  Eu 
0 E\ j)  show  a  dramatic  increase  and  follows  the  trend 

Ni/SWCNT(5,0)  <  Ni/SWCNT(10,0) 

<  Ni/SWCNT(15,0)  <  Ni/MWCNT. 

These  increases  depend  on  the  CNT  volume  fractions  of  the 
composites  and  approximately  proportional  to  the  CNT  vol¬ 
ume  fractions.  Similar  trend  is  seen  for  E22  (E22)  and  E33 
(Ej3).  For  SWCNTs,  a  slight  but  gradual  decrease  is  seen  for 
the  E  values  depending  on  the  CNT  diameter  and  volume 
fractions.  Interestingly,  E22  and  EL  increases  for  Ni/ 
MWCNT.  For  example,  in  this  case,  eJ2  (or  Ej3)  is  increased 
by  2  GPa  for  low  CNT  volume  fraction  and  an  average  of  8 
GPa  for  high  CNT  volume  fraction,  respectively,  which  sug¬ 
gests  that  Ni/CNT  composites  with  significantly  larger 
MWCNT  diameter  and  with  high  CNT  volume  fractions 
might  increase  the  Young’s  modulus  in  both  longitudinal  and 
transverse  directions.  However,  verification  of  this  hypoth¬ 
esis  is  beyond  the  scope  of  the  present  paper. 

Figure  4  shows  the  calculated  Poisson’s  ratio,  v,  for  Ni 
and  the  Ni/CNT  composites.  Our  predicted  value  for  pure  fee 
Ni  is  0.374.  Several  unique  Poisson’s  ratios  can  be  defined 
for  the  composites  depending  on  the  orthogonal  directions  of 
the  unit  vector  along  the  longitudinal  direction  and  the  trans¬ 
verse  directions  of  the  composites:  (1)  v\2  and  i^3  involve 
the  strains  between  the  longitudinal  direction  and  the  trans¬ 
verse  directions.  Here,  a  uniaxial  strain  is  applied  in  the  lon¬ 
gitudinal  direction  while  the  transverse  directions  respond  to 
this  freely.  Usually  this  indicates  contraction  in  the  trans¬ 
verse  direction  during  longitudinal  elongation;  (2)  v\2  and 
vj3  are  similar  to  v\2  and  v\3,  however,  in  this  case  a  uniaxial 
strain  is  applied  in  the  transverse  direction  while  response 
(contraction)  is  expected  in  the  longitudinal  direction;  and, 
(3)  in  4,  both  transverse  directions  are  involved,  however,  a 
uniaxial  strain  is  applied  in  one  transverse  direction  and  con¬ 
traction  is  seen  in  the  other  transverse  direction.  Note,  since 
similar  values  are  found  for  and  r^3,  as  well  as  v[2  and 
vj3,  we  show  an  average  values  for  these. 


From  Fig.  4,  the  following  trends  are  seen:  (1)  i^2  or  ^13 
show  slight  variation  across  the  different  Ni/CNTs  systems 
starting  from  Ni/SWCNT(5,0)  in  the  left  to  Ni/ 
SWCNT(15,0)  and  Ni/MWCNT  on  the  right.  In  the  compos¬ 
ites  with  high  CNT  volume  fraction,  slight  decrease  is  seen 
from  left  to  right  depending  on  CNT  diameters.  The  largest 
decrease  is  seen  for  Ni/MWCNT.  In  this  case,  for  high  CNT 
volume  fractions,  or  ^13  are  found  to  be  0.34,  a  value 
only  marginally  smaller  than  that  of  0.37  of  Ni.  (2)  In  v[2  or 
v[3,  a  significant  and  rapid  decrease  is  seen  across  for  the 
Ni/CNT  systems  studied  (left  to  right  in  Fig.  4).  Composites 
with  high  volume  fraction  show  a  far  larger  decrease.  The 
Ni/MWCNT  nanocomposites  show  larger  decrease  than  any 
of  the  Ni/SWCNT  composites.  For  example,  for  Ni/CNTs 
with  a  high  CNT  volume  fraction,  the  v[2  (or  vj3)  is  found  to 
be  0.31,  0.26,  0.22,  and  0.18,  for  Ni/SWCNT(5,0),  Ni/ 
SWCNT(10,0),  Ni/SWCNT(15,0),  and  Ni/MWCNT,  respec¬ 
tively.  An  interesting  difference  in  trend  is  seen  between  the 
longitudinal  i^2  (or  v\3)  and  the  transverse  v[2  (or  vj3).  Here, 
CNT  volume  fractions  dependency  is  seen  for  both  of  these, 
a  much  steeper  decrease  is  seen  for  v[2  (or  vj3)  across  the 
Ni/CNTs.  (3)  For  ^23,  a  large  but  gradual  increase  is  seen 
across  the  Ni/CNTs,  which  depend  on  the  CNT  diameters. 
This  increase  is  also  found  to  be  CNT  volume  fraction  de¬ 
pendent,  with  a  larger  increase  for  composites  with  high  vol¬ 
ume  fractions.  The  variations  in  v23  indicate  how  Ni/CNTs 
with  different  kind  of  CNTs  would  respond  when  an  uniaxial 
strain  is  applied  along  one  transverse  direction  while  other  is 
being  freely  contracted.  The  significantly  larger  changes  in 
v[2  (or  vj3)  and  v23  indicate  that  these  composites  have  mark¬ 
edly  different  behavior  along  longitudinal  and  transverse  di¬ 
rections  when  uniaxial  strain  is  applied. 

For  all  composites,  the  calculated  shear  modulus  (G)  for 
single  crystals  decrease  compared  to  nickel.  For  composites 
with  high  CNT  volume  fractions,  a  larger  decrease  is  seen  as 
compared  to  those  with  low  volume  fractions.  For  Ni/ 
SWCNT,  both  G23  and  G13  or  G12  decrease  with  increasing 
CNT  diameters. 

2.  Poly  crystalline  elastic  moduli 

Real  materials  are  usually  polycrystalline  aggregates  of 
randomly  oriented  single-crystal  grains,  each  exhibiting 
single-crystalline  elastic  properties.  These  materials  exhibit 
isotropic  properties.  The  elastic  moduli  of  polycrystalline  Ni/ 
CNT  aggregates  were  computed  by  averaging  appropriate 
single-crystal  constants  {ct])  using  the  expressions  of  Voigt50 
and  Reus s. 51  These  expressions  provide  upper  and  lower 
bounds  for  the  averaged  bulk  and  shear  modulus.  As  shown 
thermodynamically  by  Hill,52  the  average  of  the  Voigt  and 
Reuss  bounds  provide  better  estimates  of  the  elastic  moduli 
of  poly  crystalline  materials.  By  treating  the  poly  crystalline 
materials  as  aggregates  of  single  crystals  with  random  orien¬ 
tations,  the  isotropic  poly  crystalline  elastic  moduli  can  be 
computed  as  averages  of  anisotropic  single-crystal  elastic 
constants.  The  theoretical  lower  and  upper  bounds  to  the  true 
polycrystalline  bulk  modulus  B  and  shear  modulus  G  are 
given  by  Reuss  (assuming  uniform  stress  throughout  a  poly¬ 
crystal)  and  Voigt  (assuming  uniform  strain)  as  follows: 
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TABLE  IV.  Poly  crystalline  bulk  (B)  and  shear  (G)  moduli  calculated  (in  GPa)  for  Ni  and  nickel/CNT  composites.  The  anisotropies  of  the 
bulk  modulus  along  the  a  axis  with  respect  to  the  b  and  c  axes  are  also  given.  The  values  are  shown  for  Voigt  (V),  Reuss  (R),  and  Hill  (H) 
approximations.  Low  and  High  are  the  composites  low  and  high  CNT  volume  fractions. 


Ni(fcc) 

Gv  (95.2) 

By 

(180.7) 

Gr 

(78.04) 

Br 

(180.7) 

Gh 

(86.6) 

Bh 

(180.7) 

Low 

High 

Low 

High 

Low 

High 

Low 

High 

Low 

High 

Low 

High 

Ni/SWCNT(5,0) 

93.3 

89.6 

180.9 

181.4 

77.6 

77.0 

181.1 

181.2 

85.5 

83.3 

181.1 

181.3 

Ni/SWCNT(10,0) 

93.0 

87.7 

180.6 

180.4 

78.3 

77.5 

180.5 

179.7 

85.6 

82.6 

180.5 

180.0 

Ni/SWCNT(15,0) 

92.0 

83.7 

179.1 

174.4 

78.1 

74.4 

179.0 

173.7 

85.0 

79.1 

179.0 

174.1 

Ni/MWCNT 

93.2 

92.6 

182.2 

186.7 

80.7 

83.8 

182.0 

183.2 

86.9 

88.2 

182.0 

184.9 

Br  =  1/[SU  +  £33  +  S33  +  2(su  +  s23  +  s13)], 

Bv=  1/9  (cn  +  c33  +  c33)  +  2/9(c12  +  c23  +  c13), 

Gr  =  lS/^^n  +  S33  +  S33)  -  4 (S12  +  s23  +  ^13) 

+  3(544  +  ^55  +  ^66)1  ’ 

GV=  {C\\  +  c22  +  c33  ~  c12_  c23  ~  +  (c44  +  C55  +  C66)/ 5, 

where  are  elastic  compliances  and  their  values  can  be 
obtained  through  an  inversion  of  the  elastic  constant  matrix, 
S=C~l.  We  estimated  the  true  poly  crystalline  values  by  using 
the  Voigt  and  Reuss,  and  Hill  (VRH)  approximations:52  BH 
=  (Br+Bv)/2  and  GH=(GR+Gv)/2.  The  average  Young’s 
modulus  and  Poisson’s  ratio  for  polycrystalline  aggregates 
can  then  be  calculated  as  EpolyCYystaninQ=(9BHGH)/ (3BH 
+  Gh)  and  ^po]ycrystailine  =  (3Bh-2Gh)/ {<SBh+  2Gh). 

Table  IV  shows  the  calculated  bulk  ( BH )  and  shear  ( GH ) 
moduli  for  the  poly  crystalline  aggregates  of  Ni  and  the  Ni/ 
CNT  composites  using  VRH  approximations.  We  discuss  the 
results  obtained  using  the  Hill  average  only.  The  calculated 
Bh  for  poly  crystalline  nickel  is  identical  to  that  of  the  corre¬ 
sponding  single  crystal  and  also  agrees  well  with  experiment. 
The  bulk  modulus  of  Ni  determined  from  the  Birch- 
Murnaghan  equation  of  state  is  also  found  to  be  identical. 
For  Ni/CNT  composites,  no  significant  changes  are  seen  in 
Bh  except  for  Ni/CNT,  which  show  a  slight  increase  of  1 
GPa  and  4  GPa,  respectively,  for  low  and  high  CNT  volume 
fractions.  The  shear  modulus  which  measures  the  ability 
of  a  material  to  resist  transverse  deformations,  decreases 
slightly  for  Ni/SWCNTs  as  compared  to  pure  Ni.  For  fee 
nickel,  the  calculated  GH  is  found  to  be  87  GPa,  which 
agrees  well  with  experiment.49  For  Ni/MWCNT,  these  values 
increase  slightly  and  depend  on  CNT  volume  fractions. 

Pugh53  introduced  the  ratio  between  the  bulk  and  shear 
modulus,  Bh/Gh ,  for  poly  crystalline  phases  as  a  measure  of 
fracture/toughness  in  metals  by  considering  that  the  shear 
modulus  Gh  represents  the  resistance  to  plastic  deformation 
while  the  bulk  modulus  BH  represents  the  resistance  to  frac¬ 
ture.  A  high  Bh/Gh  value  is  associated  with  ductility  while  a 
low  value  is  associated  with  brittleness.  The  critical  value 
that  separates  ductile  and  brittle  materials  is  about  1.75.  Our 
calculated  BH/GH  for  pure  nickel  is  2.1.  For  Ni/CNT  com¬ 
posites,  similar  BH/GH  values  are  found.  By  considering  the 
bulk  modulus  as  a  measure  of  the  average  bond  strength  and 


the  shear  modulus  as  a  measure  of  the  resistance  to  a  change 
in  bond  angle  by  an  external  force,  Tanaka  et  al.5A  proposed 
that  Gh/Bh  represents  the  relative  directionality  of  the  bond¬ 
ing  in  the  material.  It  can  be  seen  that  there  is  a  slight  de¬ 
crease  in  these  ratios  for  Ni/SWCNTs.  Composites  with  high 
CNT  volume  fractions  show  a  greater  decrease.  For  Ni/ 
MWCNT,  there  is  no  observable  change  with  pure  Ni.  The 
Bh/Gh  and  the  GH/BH  values  indicate  that  CNT  do  not  sig¬ 
nificantly  alter  the  above-mentioned  properties  as  compared 
to  fee  Ni. 

The  average  Young’s  modulus  for  poly  crystalline  aggre¬ 
gates  (£p0iycrystaiiine)  calculated  from  the  Hill’s  formula  for 
pure  Ni  and  the  Ni/CNTs  are  presented  graphically  in  Fig.  3. 
The  Young’s  modulus  in  longitudinal  E00l  ( En )  and  trans¬ 
verse  directions,  E0l0,El00  (E22,E33)  determined  for  single¬ 
crystal  specimens  are  also  compared.  For,  Ni  the  value  cal¬ 
culated  using  Hill’s  average  is  222  GPa.  This  value  is  in 
excellent  agreement  with  the  experimental  value  of 
213.1  ±21.2  GPa,  obtained  for  coarse-grained  poly  crystal¬ 
line  aggregates.49  In  Ni/SWCNT  composites,  £poiycrystaiiine 


FIG.  5.  (Color  online)  Young’s  modulus  in  longitudinal  direc¬ 
tion,  En ,  (in  GPa)  simulated  and  calculated  using  the  rule-of- 
mixture  (RoM)  for  Ni/SWCNT  and  Ni/MWCNT  composites  with 
low-  and  high  CNT  volume  fractions.  The  straight  lines  for  each 
CNT  type  are  obtained  by  fitting  the  RoM  data  at  low  and  high 
CNT  volume  fractions.  The  simulated  values  for  each  composite  are 
presented  with  symbols. 
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TABLE  V.  The  isotropic  bulk  modulus  (B  in  GPa)  and  directional  bulk  modulus  along  the  orthorhombic 
crystallographic  axes  a,  b ,  and  c  ( Ba ,  Bb,  and  Bc )  for  pure  Nickel  (fee)  and  nickel/CNT  composites.  Low  and 
High  are  the  composites  low  and  high  CNT  volume  fractions. 


Ni(fcc) 

B  (180.7) 

Ba 

(542.1) 

Bb 

(542.1) 

Bc 

(542.1) 

Low 

High 

Low 

High 

Low 

High 

Low 

High 

Ni/SWCNT(5,0) 

181.1 

181.2 

560.8 

599.8 

530.9 

516.3 

538.8 

522.5 

Ni/SWCNT(10,0) 

180.5 

179.2 

573.3 

650.5 

527.4 

498.3 

526.6 

495.1 

Ni/SWCNT(15,0) 

179.0 

173.7 

572.6 

629.5 

520.4 

479.7 

521.2 

479.9 

Ni/MWCNT 

182.0 

183.2 

629.3 

817.4 

508.6 

462.2 

514.5 

482.9 

decreases  with  increasing  CNT  diameter.  Composites  with 
higher  CNT  volume  fractions  show  larger  decreases.  Ni/ 
MWCNT  does  not  show  an  appreciable  change  in 
^poiycrystaiiine  f°r  the  composite  with  low  volume  fraction  and 
a  slight  increase  for  composites  with  high  volume  fractions 
of  CNT.  Generally,  £p0iycrystaiiine  *s  much  larger  in  magnitudes 
than  the  single  crystal  E.  From  these  trends,  we  can  conclude 
that  embedding  SWCNTs  in  the  Ni  matrix  will  decrease  the 
poly  crystalline  Young’s  modulus  Epoiycrystaiiine-  On  the  con¬ 
trary,  the  poly  crystalline  Young’s  modulus  might  increase 
when  MWCNTs  are  embedded  in  the  Ni  matrix. 

Our  calculated  Poisson’s  ratio  for  poly  crystalline  nickel  is 
0.293.  This  is  in  reasonable  agreement  with  the  experimental 
value  of  0.31. 55  Surprisingly,  we  find  similar  values  for  Ni/ 
CNT  composites.  For  example,  for  Ni/SWCNT(5,0),  Ni/ 
SWCNT(10,0),  Ni/SWCNT(15,0),  and  Ni/MWCNT  compos¬ 
ites  with  higher  CNT  fractions,  the  calculated  Poisson’s 
ratios  are  0.301,  0.301,  0.303,  and  0.317,  respectively.  They 
are  0.296,  0.295,  0.295,  and  0.294  for  the  same  composites 
with  lower  CNT  volume  fractions.  The  calculated  Poisson’s 
ratios  indicate  that  the  Ni/CNT  composites  remain  similarly 
stable  against  shear  as  compared  to  pure  Ni. 

3.  Elastic  anisotropy 

Elasticity  describes  the  response  of  a  crystal  under  exter¬ 
nal  strain  and  provides  key  information  about  the  bonding 
characteristic  between  adjacent  atomic  planes  and  the  aniso¬ 
tropic  character  of  the  solid.26  The  shear  anisotropic  factors 
measure  the  degree  of  anisotropy  in  the  bonding  between 
atoms  in  different  planes.  Anisotropy  in  orthorhombic  mate¬ 
rials  arises  from  the  anisotropy  in  their  shear  and  linear  bulk 
moduli. 

The  shear  anisotropic  factor  is  defined  as 


Aj  -  4c44/(cn  +  c33  -  2 cl3) 

for  the  {100}  shear  planes  in  (Oil)  and  (010)  directions, 


A2  -  4c55/(c22  +  c33  -  2c23) 

for  the  {010}  shear  planes  in  (101)  and  (001)  directions,  and 


A3  -  4c66/(cn  +  c22  -  2c if) 

for  the  {001}  shear  planes  in  (110)  and  (010)  directions.  For 
an  isotropic  crystal,  the  factors  A1?  A2,  and  A3  must  be  one. 
The  deviation  from  one  is  a  measure  of  the  degree  of  the 
elastic  anisotropy  of  crystals. 


The  anisotropy  of  the  bulk  modulus  along  the  a  axis  with 
respect  to  b  and  c  axes  can  be  estimated  using  the  following 
equations: 


Ba 

Bh 


Ba 

Br 


Note  that  a  value  of  unity  for  these  two  fractions  indicates 
elastic  isotropy.  In  the  above  equations,  Ba ,  Bb ,  and  Bc  are 
the  bulk  moduli  along  different  crystal  axes,  are  defined  as 


Bn  -  a 


dP 
da  ’ 


Bh  =  b 


dP 
db 9 


and 


Br  = 


dP 
c — . 
dc 


A  practical  measure  of  elastic  anisotropy  for  poly  crystalline 
materials  is  defined  as  follows: 


(By-BR) 

(By  +  BR) 


and 


(Gy~  GR ) 
(Gy+  GR ) 


where  B  and  G  are  the  bulk  and  shear  moduli  and  the  sub¬ 
scripts  V  and  R  denote  the  Voigt  and  Reuss  limits.  For  these 
expressions,  a  value  of  zero  identifies  elastic  isotropy  and  a 
value  of  1  (100%)  is  the  largest  possible  anisotropy. 

The  bulk  modulus  and  the  directional  bulk  modulus  are 
given  in  Table  V.  Note  that  the  CNTs  are  aligned  along  the  a 
axis.  It  is  clear  that  Ba  is  higher  than  Bb  and  Bc ,  which  are 
nearly  identical.  In  the  composites,  Ba  increases  with  volume 
fraction  of  CNT  while  Bb  and  Bc  both  decrease  as  compared 
to  Ni.  The  directional  bulk  modulus  of  Ni  is  542  GPa  along 
all  three  axes.  For  Ni/MWCNT  composites,  the  highest  (817 
GPa)  and  lowest  (462  GPa)  values  for  Ba  and  Bb  or  Bc  are 
seen  for  high  CNT  volume  fractions.  For  Ni/SWCNTs,  de¬ 
pendencies  on  the  diameters  are  also  seen. 

The  shear  anisotropic  factors  are  given  in  Table  VI.  For 
fee  Ni,  we  calculated  a  shear  anisotropic  factor  As  is  2.52 
while  the  value  reported  in  the  literature56  value  is  2.63.  This 
indicates  significant  deviation  from  isotropy.  We  also  see  that 
the  shear  anisotropic  factors  A1?  A2,  and  A3  for  the  compos¬ 
ites  depend  on  CNT  concentrations.  Similar  changes  are  seen 
for  A 1  and  A3  where  these  values  are  found  to  be  lower  than 
pure  nickel.  A  steeper  decrease  for  these  is  seen  for  compos¬ 
ite  with  high  CNT  volume  fractions.  For  SWCNTs,  a  larger 
decrease  is  seen  for  CNTs  with  larger  diameter.  A2  for  the 
composites  remains  similar  to  that  for  pure  nickel. 
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TABLE  VI.  The  shear  anisotropic  factors  A  1?  A2,  A3,  and  AG  (in  %),  AB  (in  %)  for  pure  Nickel  (fee)  and  nickel/CNT  composites  with  low 
and  high  CNT  volume  fractions. 


Ni(fcc) 

A\  (2.516) 

A2  (2.516) 

A3  (2.516) 

Ab 

(0.00) 

Ag  (9.875) 

(i) 

(i) 

Low 

High 

Low 

High 

Low 

High 

Low 

High 

Low 

High 

Low 

High 

Low 

High 

Ni/SWCNT(5,0) 

2.452 

2.155 

2.485 

2.409 

2.370 

2.026 

0.005 

0.047 

9.237 

7.542 

1.056 

1.162 

1.040 

1.148 

Ni/SWCNT(10,0) 

2.247 

1.630 

2.504 

2.452 

2.309 

1.838 

0.015 

0.180 

8.602 

6.158 

1.087 

1.305 

1.088 

1.314 

Ni/SWCNT(15,0) 

2.113 

1.269 

2.525 

2.539 

2.228 

1.657 

0.020 

0.209 

8.145 

5.874 

1.101 

1.315 

1.098 

1.311 

Ni/MWCNT 

2.002 

1.304 

2.444 

2.181 

2.000 

1.254 

0.102 

0.929 

7.218 

4.985 

1.237 

1.768 

1.223 

1.692 

From  the  anisotropy  of  the  directional  linear  bulk  modu¬ 
lus  it  is  clear  that  fee  nickel  does  not  show  any  directionality. 
However,  due  to  the  alignment  of  the  CNTs,  the  nanocom¬ 
posites  show  differences  in  the  bulk  moduli  along  different 
axes.  For  example,  along  the  direction  of  CNT  alignment 
(i.e.,  a  axis),  there  is  an  increase  in  the  bulk  modulus 
whereas  in  the  transverse  directions  (i.e.,  b  and  c  axes)  de¬ 
creases  are  seen  compared  to  Ni.  These  changes  are  more 
prominent  at  higher  CNT  concentrations.  These  results  imply 
that  the  composites  are  less  compressible  along  the  direction 
of  the  CNT  axis  and  more  compressible  along  the  transverse 
directions.  The  anisotropy  of  the  bulk  modulus  ABb  and  AB ^ 
along  the  a  axis  with  respect  to  b  and  c  axes  shows  that 
anisotropy  increase  with  the  increasing  CNT  diameter  and 
CNT  volume  fraction.  For  Ni,  the  bulk  anisotropy  (Afi)  is 
zero  indicating  no  anisotropic  compressibility.  For  compos¬ 
ites,  Ab  increases  with  increase  in  the  CNT  volume  fraction 
and  CNT  diameters  indicating  an  increase  in  compressibility. 
The  maximum  value  for  AB  is  seen  for  the  Ni/MWCNT  com¬ 
posite. 

Ag  gives  the  relative  magnitude  of  the  actual  elastic  an¬ 
isotropy  possessed  by  a  crystal.56  In  the  case  of  shear  aniso¬ 
tropy  (Ag),  ME  AM  Ni  shows  the  maximum  value  of  9.88 
while  a  value  of  9.76  is  reported  in  the  literature.56  This 
suggests  that  Ni  has  large  shear  anisotropy.  A  decrease  is 
seen  for  composites  with  larger  CNT  radii  and  volume  frac¬ 
tions.  Thus,  it  can  be  concluded  that  the  composites  are  also 
elastically  anisotropic.  However,  there  are  significant  effects 
of  CNTs  within  the  composites,  the  AG  decrease  and  larger 
decreases  are  seen  for  composites  with  higher  CNT  volume 
fractions. 

It  is  interesting  to  note  that  the  composites  have  the  high¬ 
est  directional  bulk  modulus  along  the  direction  in  which  the 
nanotubes  are  aligned.  In  the  transverse  directions,  the  values 
are  comparable.  This  indicates  that  the  compressibility  is 
smallest  along  the  longitudinal  direction  and  largest  along 
the  transverse  direction.  Along  the  CNT  alignment  direction, 
the  bulk  modulus  increases  compared  to  pure  nickel  and  it 
decreases  in  the  transverse  directions.  CNT  concentration 
also  affects  the  bulk  modulus  significantly. 

4.  Interfacial  characteristics 

We  have  closely  monitored  the  motion  of  nickel  and  car¬ 
bon  atoms  in  Ni/SWCNT(10,0)  with  low  CNT  volume  frac¬ 
tion  at  0.0%  and  0.2%  applied  strains.  We  find  both  nickel 
and  carbon  atoms  move  in  concert  with  each  other  in  the 


direction  of  applied  strain.  There  is  little  or  no  movement  of 
atoms  in  the  transverse  directions,  except  for  a  few  atoms  in 
the  central  portion  of  the  cell  where  Ni  atoms  are  found  to 
perturb  laterally  in  the  perpendicular  directions.  In  the  vicin¬ 
ity  of  the  nanotube,  it  is  seen  that  first  few  Ni  atoms  move 
away  and  others  shift  toward  the  center  during  minimization 
in  this  region.  The  overall  relative  shifts  of  atomic  positions 
during  the  loading  can  be  quantitatively  estimated  using  the 
root  mean-square  deviations  (RMSDs)  of  the  atoms  in  the 
simulation  cell  from  those  two  strain  limits.  It  is  seen  that  the 
RMSD  for  Ni/CNT  system  is  0.0214  A  in  the  longitudinal 
direction,  which  interestingly,  is  virtually  identical  to  the 
RMSDs  calculated  for  CNT  and  Ni  atoms  separately.  Note 
that  these  RMSD  calculations  were  performed  by  taking  on 
the  central  portion  of  the  original  cells  but  reducing  it  by  \ 
X  ^  X  2  in  x,y,  and  z  directions.  In  the  transverse  directions, 
the  RMSD  is  0.0007  A,  indicating  very  small  or  no  relative 
movement  due  to  applied  strain.  These  results  suggest  that 
both  carbon  and  nickel  atom  move  in  concert  during  uniaxial 
loading  in  the  longitudinal  direction,  without  sliding  of  C 
atoms  relative  to  the  nickel  atoms.  This  is  possible  in  the 
case  of  a  relatively  strong  Ni/CNT  interface.  This  finding  is 
also  supported  by  our  analysis  of  rule-of-mixture  (RoM)  re¬ 
sults,  discussed  in  the  next  section. 

5.  Analysis  of  rule-of-mixturc 

The  properties  of  fiber-reinforced  composites  have  been 
extensively  studied  to  predict  their  elastic  properties.57  One 
such  approach  is  the  RoM.58,59  For  a  CNT-reinforced  nano¬ 
composites,  under  uniaxial  loading,  the  dependence  of  the 
Young’s  modulus  on  the  CNT  volume  fraction  can  be  theo¬ 
retically  estimated  as 

Ec=VfEf+(l-Vf)Em, 

where  Ec  is  the  predicted  Young’s  modulus  of  the  composite, 
Ej-  is  the  Young’s  modulus  of  the  fiber  reinforcement  (CNT), 
Em  the  Young’s  modulus  of  the  matrix,  and  Vj-  is  the  volume 
fraction  of  the  CNTs. 

We  have  estimated  the  Young’s  modulus  in  the  longitudi¬ 
nal  direction  (. En )  for  all  Ni/CNT  composites  and  compared 
them  with  the  results  obtained  directly  from  simulation.  Fig¬ 
ure  5  shows  En  as  a  function  of  Vj  for  Ni/CNT  composites 
(here  the  full  lines  plots  of  the  standard  RoM  expression).  It 
is  clear  that  the  RoMs  fit  with  the  simulated  results  well  and 
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implicates  a  strong  Ni/CNT  interfacial  interaction.  The  RoM 
En  increases  linearly  with  volume  fractions  and  also  shows 
dependence  with  CNT’s  radii. 

Comparison  of  the  RoM  estimated  En  with  simulated  val¬ 
ues  in  the  case  of  SWCNTs  shows  an  opposite  trend.  Accord¬ 
ing  to  the  RoM,  we  should  expect  SWCNT(15,0)  to  give  the 
largest  value  of  En  compared  to  that  of  SWCNT(5,0).  How¬ 
ever,  simulated  results  predict  the  opposite.  SWCNT(10,0) 
shows  a  more  remarkable  agreement  between  calculated  and 
simulated  values  than  any  other  CNTs.  Note  that  the  plots  in 
Fig.  5  are  very  useful  as  these  can  be  used  to  predict  the 
modulus  for  CNT  with  any  diameter  and  volume  fraction. 

To  reasonably  compare  the  different  CNT  composites,  one 
can  use  the  rate  of  increase  in  Young’s  modulus  with  volume 
fraction  ( dEn/dVf )  as  a  convenient  yardstick  for  assessing 
the  effect  of  reinforcement  (CNT  effectiveness).  This  useful 
parameter  shows  the  magnitude  of  the  increase  in  stiffness 
and  the  CNT  content  required  achieving  that  increase.  For  all 
four  composites,  the  simulated  En  increases  approximately 
linearly  with  volume  fraction  of  CNT.  The  dEn/dVf  are  707 
GPa,  629  GPa,  474  GPa,  and  951  GPa  for  Ni/SWCNT(5,0), 
Ni/SWCNT(10,0),  Ni/SWCNT(15,0),  and  Ni/MWCNT,  re¬ 
spectively.  The  highest  value  of  dEn/dVf  (951  GPa)  was 
found  with  the  MWCNT.  Using  the  RoM  data,  the  enhance¬ 
ment  dEn/dVf  is  found  to  be  415  GPa,  695  GPa,  755  GPa, 
and  603  GPa,  respectively,  for  Ni/SWCNT(5,0),  Ni/ 
SWCNT(10,0),  Ni/SWCNT(15,0),  and  Ni/MWCNT.  This 
opposite  trend  between  the  RoM  and  simulated  results  sug¬ 
gests  that  the  reinforcement  mechanism  is  not  only  from 
mixing  alone.  This  deviation  likely  arises  from  the  large 
surface-area-to-volume  ratios  of  CNTs  that  affect  matrix 
properties  beyond  what  would  be  predicted  from  RoM  cal¬ 
culations. 

It  is  known  that  the  performance  of  a  reinforced  compos¬ 
ite  depends  on  the  interfacial  characteristics  of  the  reinforce¬ 
ment  fillers  and  the  matrix  materials.60  For  example,  strong 
adhesion  forces  at  the  interface  lead  to  superior  mechanical 
properties.  It  also  critically  depends  on  the  effectiveness  of 
the  interfacial  stress  transfer,  which  in  turn  depends  on  the 
nature  and  strength  of  the  interface.  Chemical  bonding,  na¬ 
nomechanical  interlocking,  and  nonbonded  (van  der  Waals 
and  electrostatic)  interactions  yield  strong  interaction  be¬ 
tween  a  polymer  matrix  and  CNTs.  Furthermore,  chemical 
functionalization  of  the  filler  surface  and  increasing  non¬ 
bonded  attractive  interactions  between  the  filler  and  the  ma¬ 
trix  increases  such  interfacial  adhesion  and  improve  the  me¬ 
chanical  properties  of  composites.  For  most  polymers,  the 
polymer/CNT  interfacial  stress  transfer  is  low  with  values 
typically  less  than  50  MPa.61 

For  metal/CNT  interfaces,  above-discussed  information  is 
not  available.  The  large  increase  in  En  for  Ni/CNT  compos¬ 
ites  indicates  that  there  exists  a  surprisingly  strong  adhesive 
force  between  the  CNT  and  the  Ni  matrix.  An  applied  load  is 
largely  transferred  to  the  reinforcing  CNTs,  indicating  an  ef¬ 
fective  stress  transfer  across  the  Ni-CNT  interfaces  that  in 
turn  enhance  the  mechanical  response.  We  would  like  to  em¬ 
phasize  that  a  thorough  investigation  using  high-level  theo¬ 
ries  is  necessary  to  fully  understand  the  chemistry  at  the 
Ni-CNT  interfaces.  Recent  first-principles  calculations  on 
different  structural  models  (e.g.,  top  fee,  bridge  top,  etc.)  of 


graphene  on  Ni(lll)  show  very  weak  or  no  binding  interac¬ 
tion  between  graphene  and  Ni  surface,62,63  which  may  be  due 
to  shortcomings  of  current  levels  of  density-functional  theory 
utilized  in  solid-state  calculations.  Simulations  using  novel 
and  adequate  theory  to  rationalize  the  interfacial  behavior  in 
Ni/CNT  composites  are  beyond  reach  at  this  time  due  to  the 
sheer  size  of  these  systems. 

IV.  CONCLUSIONS 

A  MEAM  formalism  has  been  used  to  investigate  the  me¬ 
chanical  properties  of  Ni/CNT  nanocomposites  made  by  in¬ 
corporating  SWCNTs  with  three  different  diameters  and  a 
MWCNT  into  fee  Ni  matrices.  Two  different  volume  frac¬ 
tions  for  each  CNT  were  used  to  assess  the  CNT  concentra¬ 
tion  effects  upon  mechanical  properties  of  the  composites. 
Using  quasistatic  energy-minimization  methods  and  periodic 
boundary  conditions,  the  single-crystal  elastic  constants, 
Young’s  and  shear  moduli,  Poisson’s  ratios  for  pure  nickel 
matrix,  and  Ni/CNT  composites  were  calculated.  The  elastic 
anisotropy  in  the  composites  was  also  systematically  inves¬ 
tigated.  The  volume  fraction  dependence  of  the  Young’s 
modulus  was  analyzed  using  rule-of-mixture  and  compared 
that  with  the  simulated  values. 

Calculations  show  that  pristine  SWCNTs  have  high 
Young’s  modulus  values  of  0.55  TPa,  0.83  TPa,  and  0.89  TPa 
for  the  SWCNT(5,0),  SWCNT(10,0)  and  SWCNT(15,0),  re¬ 
spectively.  The  MWCNT  made  from  the  above  three 
SWCNTs  show  a  Young’s  modulus  of  0.74  TPa.  These  val¬ 
ues  are  in  reasonable  agreement  with  the  observed  and  theo¬ 
retically  predicted  values  published  in  the  literature. 

The  results  for  Ni/CNT  composites  are  primarily  com¬ 
pared  with  the  results  for  a  pure  fee  nickel  matrix.  Our  cal¬ 
culations  predict  that  all  Ni/CNTs  satisfy  the  conditions  for 
mechanical  stability.  In  the  cases  of  single-crystal  Ni/CNT 
composites,  the  Young’s  modulus  shows  a  significant  in¬ 
crease  in  the  longitudinal  direction  (the  direction  of  CNT 
alignment)  and  a  moderate  decrease  in  the  transverse  direc¬ 
tions.  It  is  found  that  CNT  volume  fractions  have  large  ef¬ 
fects  on  the  single-crystal  Young’s  modulus.  For  pure  fee 
nickel,  our  calculated  Young’s  modulus  is  137  GPa.  In  the 
composites  with  high  CNT  volume  fraction,  the  single¬ 
crystal  Young’s  modulus  along  the  longitudinal  direction  is 
found  to  be  160  GPa,  183  GPa,  and  201  GPa  for  Ni/ 
SWCNT(5,0),  Ni/SWCNT(10,0),  and  Ni/SWNCT(15,0),  re¬ 
spectively.  Ni/MWCNT  shows  the  highest  value  of  266  GPa 
for  Young’s  modulus  in  same  direction  when  high  CNT  vol¬ 
ume  fractions  are  used.  In  the  transverse  directions,  the 
Young’s  modulus  decreases  by  1-7  GPa  for  Ni/SWCNTs  de¬ 
pending  on  the  CNT  volume  fraction  and.  For  Ni/MWCNT 
however,  3-12  GPa  increase  is  seen. 

For  the  composites  with  SWCNTs,  the  poly  crystalline 
Young’s  modulus  decreases  with  increasing  CNT  diameters 
and  CNT  volume  fractions.  In  the  Ni/MWCNT  nanocompos¬ 
ite,  however,  a  slight  increase  is  seen  at  high  CNT  volume 
fractions,  indicating  possibility  of  fabricating  composites 
with  enhanced  overall  mechanical  properties  when 
MWCNTs  with  sufficiently  larger  size  and  volume  fractions 
were  used.  The  polycrystalline  shear  modulus  remain  similar 
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to  pure  nickel  for  the  composites  with  Ni/SWCNTs,  how¬ 
ever,  slight  increase  is  seen  for  Ni/MWCNT  composite  with 
higher  CNT  volume  fractions.  The  result  indicates  that  Ni/ 
CNT  composites  with  embedded  single-walled  nanotubes 
with  large  CNT  volume  fraction  would  give  slightly  inferior 
mechanical  properties  than  that  of  pure  nickel.  However, 
slightly  enhanced  mechanical  properties  may  be  expected 
when  MWCNT  with  sufficiently  large  diameter  and  volume 
fractions  are  used. 

Due  to  aligned  CNTs,  the  composites  show  difference  in 
the  bulk  modulus  along  different  axes.  The  present  calcula¬ 
tions  show  that  the  bulk  modulus  along  the  longitudinal  di¬ 
rection  is  much  higher  than  the  transverse  directions  and  in¬ 
creases  with  the  CNT  volume  fractions.  The  shear 
anisotropic  factors  (A ,,  A2,  and  A3 )  for  the  composites  also 
depend  on  the  CNT  concentrations  and  are  found  to  depend 
on  the  directions  of  CNT  alignment.  The  bulk  anisotropy 
(Afi)  increases  with  increase  in  the  CNT  volume  fraction  and 
CNT  diameter,  indicating  an  increase  in  compressibility.  The 
maximum  value  for  AB  is  seen  for  Ni/MWCNT  nanocompos¬ 
ite.  A  decrease  is  seen  in  the  shear  anisotropy  for  composites 
with  larger  CNT  radii  and  volume  fractions. 
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The  observations  that  (1)  the  Young’s  moduli  obtained 
from  the  atomistic  simulations  and  rule-of-mixture  agree 
very  well,  (2)  the  reinforcement  values,  estimated  from  the 
rate  of  increase  in  the  longitudinal  Young’s  modulus  as  a 
function  of  CNT  volume  fraction,  are  remarkably  high,  and, 
(3)  the  Ni  and  C  atoms  at  the  Ni/CNT  interface  move  in 
concert,  in  the  direction  of  applied  strain,  suggest  a  very 
strong  Ni/CNT  interface  in  the  Ni/CNT  composites. 
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